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ABSTRACT 

The X-ray light-curves of 9 Swift XRT afterglows (050126, 050128, 050219A, 050315, 050318, 050319, 
050401, 050408, 050505) display a complex behaviour: a steep t~ 3 0±0 - 3 decay until ~ 400 s, followed 
by a significantly slower t~ 0-65±0.20 f a ii_ g ) which at 0.2—2 d after the burst evolves into a !-7±0.5 
decay. We consider three possible models for the geometry of relativistic blast-waves (spherical out- 
flows, non-spreading jets, and spreading jets), two possible dynamical regimes for the forward shock 
(adiabatic and fully radiative), and we take into account a possible angular structure of the outflow 
and delayed energy injection in the blast-wave, to identify the models which reconcile the X-ray 
light-curve decay with the slope of the X-ray continuum for each of the above three afterglow phases. 
By piecing together the various models for each phase in a way that makes physical sense, we identify 
possible models for the entire X-ray afterglow. The major conclusion of this work is that a long-lived 
episode of energy injection in the blast-wave, during which the shock energy increases at t 1 ' ^ ' 5 , 
is required for 5 afterglows and could be at work in the other 4 as well. For some afterglows, there 
may be other mechanisms that can explain the t < 400 s fast falling-off X-ray light-curve (e.g. the 
large-angle GRB emission), the 400 s— 5 h slow decay (e.g. a structured outflow), or the steepening at 
0.2— 2d (e.g. a jet-break, a collimated outflow transiting from a wind with a r~ 3 radial density profile 
to a homogeneous or outward-increasing density region). Optical observations in conjunction with the 
X-ray can distinguish among these various models. Our simple tests allow the determination of the 
location of the cooling frequency relative to the X-ray domain and, thus, of the index of the electron 
power-law distribution with energy in the blast-wave. The resulting indices are clearly inconsistent 
with an universal value. 

Key words: gamma-rays: bursts - ISM: jets and outflows - radiation mechanisms: 
non-thermal - shock waves 



1 INTRODUCTION 

Pre-Swift observations of Gamma-Ray Burst (GRB) after- 
glows have led to great strides in their theoretical interpre- 
tation, while leaving some major unanswered questions. 

The radio, optical, and X-ray emission of GRB after- 
glows exhibit a power-law decrease with time (F v oc t~ a ) 
from hours to tens of days after the burst, with the tem- 
poral index a consistent with the slope (3 of the power-law 
continuum (F u oc v" 13 ) within the framework of relativistic 
spherical blast-waves (Meszaros & Rees 1997) or of spread- 
ing relativistic jets (Rhoads 1999). The collimation of the 
GRB ejecta yields a steepening of the power-law decay when 
the relativistic beaming has decreased sufficiently that the 
jet boundary becomes visible. Such a steepening has been 
observed for the first time in the optical light-curve of the 
afterglow 990123 (Kulkarni et al. 1999). Since then about 10 
other afterglows have displayed an optical light-curve break 
at about 1 day after the burst. The achromaticity of a jet 
light-curve break has not been clearly proven by pre-Swift 



observations because the X-ray light-curves were not mon- 
itored over a time long enough to capture the jet-break. 
Furthermore, the radio light-curves were usually poorly sam- 
pled during the first day and strongly affected by interstellar 
scintillation. The observations of many X-ray afterglows by 
Swift, together with ground-based optical observations, will 
enable us to test achromaticity of the afterglow light-curve 
break, as appears to be the case for the afterglow 050525A 
(Blustin et al. 2005). 

The quenching of the interstellar scintillation of the ra- 
dio afterglow 970508 (Frail, Waxman & Kulkarni 2000) has 
confirmed that the source size increases as expected for a rel- 
ativistic blast-wave, providing another test for this model. 
The decrease of the scintillation has also been observed in 
the radio afterglows 991208, 021004, and 030329. Further 
testing has been prompted by the detection of the optical 
afterglows of GRBs 990123 and 021211 at very early times 
(Akerlof et al. 1999, Fox et al. 2003, Li et al. 2003), starting 
at about 100 s after the burst. The steep decays (a = 1.8 
and 1.6, respectively) exhibited by these afterglows in the 
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first 20 minutes can be attributed to the GRB ejecta en- 
ergized by internal shocks (Meszaros & Rees 1997, 1999) 
or by the reverse shock*caused by the interaction with the 
circumburst medium (Sari & Piran 1999). 

A significant discrepancy between afterglow observa- 
tions and theoretical expectations exists for the radio after- 
glows of GRBs 991208, 991216, 000301c, and 010222, whose 
decay over 1-2 decades in time is substantially slower than 
that of the optical emission (Frail et al. 2004, Panaitescu & 
Kumar 2004). A change in the blast-waves dynamics, such as 
the transition to semi-relativistic dynamics, is not a possible 
explanation, because the different radio and optical decays 
are observed over time ranges which overlap substantially. 
Our analysis (Panaitescu & Kumar 2004) of these afterglows 
shows that evolving microphysical parameters cannot de- 
couple the optical and radio decays. This decoupling may 
be achieved if there is an extra radio emission arising from 
some late ejecta, energized by a reverse shock. For the optical 
afterglow to remain unaffected, the incoming ejecta should 
not alter the dynamics of the blast-wave, i.e. they should 
carry less kinetic than that already existing in the swept-up 
circumburst medium. 

The Swift measurements of the X-ray afterglow emis- 
sion, starting from 100 s after the burst, opens new pos- 
sibilities for testing the blast-wave model and for refining 
its details. The XRT 0.2-10 keV light-curves of the 9 X- 
ray afterglows (050126, 050128, 050219A, 050315, 050318, 
050319, 050401, 050408, 050505) presented by Campana et 
al. (2005), Chincarini et al. (2005), and Tagliaferri et al. 
(2005), have shown that some X-ray afterglows decay very 
fast (F x oc t~ [i ) within the first few minutes after the burst, 
as reported previously for the afterglows 990510 (Pian et 
al. 2001) and 010222 (in't Zand et al. 2001), followed by a 
slower decay phase (F x oc t~ 2 ^ 3 ), and a break to a steeper 
decay (F x oc t^ 5 '' 3 ) at a later time, ranging from 1 hour to 
1 day. The purpose of this paper is to investigate what fea- 
tures of the blast-wave model are required to accommodate 
the various decays of these Swift X-ray afterglows. 

Barthelmy et al. (2005) have shown that the very early 
fast decay of the X-ray emission of the afterglows 050315 
and 050319 can be understood as the GRB emission from 
the fluid moving at angles larger than the inverse of the 
forward shock's Lorentz factor. Due to the curvature of the 
emitting surface, this large angle emission arrives at the ob- 
server at an ever increasing time and ever decreasing fre- 
quency (Kumar & Panaitescu 2000). However, for two other 
Swift afterglows with an early, fast decaying X-ray emis- 
sion (050126 and 050219A), Tagliaferri et al. (2005) have 
found that the 15-350 keV GRB emission extrapolated to 
the XRT 0.2-10 keV band (under the assumption that the 
burst power-law spectrum extends unbroken to lower ener- 
gies) falls short of the flux measured at the beginning of the 
X-ray observations. If the burst spectrum has a break below 
15 keV, below which it is harder, then the GRB extrapolated 
flux would be even less. This suggests that the fast falling-off 



* Later it became clear that the latter interpretation is not so 
straightforward: the reverse-shock microphysical parameters re- 
quired to accommodate the early optical light-curves of the after- 
glows 990123 and 021211 imply a magnetized ejecta (Fan et al. 
2002, Kumar & Panaitescu 2003, Zhang, Kobayashi & Meszaros 
2003) 



X-ray emission does not arise from the same mechanism as 
the burst itself, and that it may arise in the forward shock. 
Therefore, for at least these two last bursts, we shall test 
whether the very early X-ray emission can have the same 
origin as the rest of the afterglow. 

The steepening observed at later times is most naturally 
attributed to a collimated outflow (jet), hence we shall test 
if the pre- and post-break X-ray light-curve indices and the 
spectral slopes are consistent with this interpretation. 



2 THE X-RAY LIGHT-CURVE DECAY INDEX 

For the dynamics and collimation of the relativistic blast- 
wave, we consider three cases: i) a spherical GRB remnant, 
in the sense that observations were done at a time when 
the afterglow Lorentz factor V was larger than the inverse 
of the jet opening 9j and, hence, the effects associated with 
collimation were not yet detectable, ii) a jet whose edge is 
visible (F6j < 1) and which does not expand laterally (be- 
cause it is embedded in an outer outflow, but whose emission 
is dimmer), and Hi) a jet with sharp edges, which spreads 
laterally and is observed when T9j < 1. These models will 
be named S , j , and J , respectively. 

At a frequency above that of the synchrotron peak, 
the index a of the light-curve power-law decay depends on 

i) the index p of the power-law electron distribution with 
energy 

dN _ v 

ii) the density stratification of the circumburst medium 
(CBM), for which we assume a power-law profile 



n(r) oc ' 



s < 3 



(2) 



which comprises a homogeneous CBM (s = 0) and a pre- 
ejected wind at constant speed and mass-loss rate (s = 2), 
the condition s < 3 being required for a decelerating blast- 
wave, 

in) the location of the cooling frequency u c relative to the 
observing band. The v c is the synchrotron characteristic fre- 
quency corresponding to an electron energy for which the ra- 
diative (synchrotron + inverse Compton) timescale is equal 
to the electron age. 

The expressions for a(p,s) for the S model are given in 
Meszaros & Rees (1997) and Sari, Narayan & Piran (1998) 
for s = and in Chevalier & Li (2000) for s = 2. Rhoads 
(1999) and Sari, Piran & Halpern (1999) have shown that, 
for the J model, a — p, irrespective of the location of v c and 
CBM stratification. These and other results for the S and 
j models are summarized below. 

2.1 Adiabatic Afterglows 

Because we will determine from observations the required 
structure of the CBM, i.e. the parameter s, we start from 
the most general expressions for the evolution with observer 
time t of the afterglow spectral properties: peak flux F p , and 
frequencies z/» and v c . As derived in Panaitescu & Kumar 
(2004), they are: 

d In i/i 3 d In v c 3s — 4 
dint _ _ 2 ' dint _ 8 - 2s 



(S,j) 



(3) 
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for both the S and j models (note that the evolution of i>i is 
in independent of s and that v c increases for s < 4/3, but 
decreases for s > 4/3), 



(S) 



dlnFp 



dint 8 -2s 

for the S model and 

dlnF D 6-s 



00 



dint 



8 -2s 



(4) 



(5) 



for the j model. For the latter, the faster decay is due to 
that the jet area is a factor (T9j) 2 smaller than that visible 
to the observer in the case of a spherical outflow. 

The synchrotron afterglow continuum is F v oc v~ fi 
(Sari, Narayan & Piran 1998), where 



r 1/2 v c < v < Vi 

(3=1 (p- l)/2 vi < v < v c 
\ p/2 Vi, v c <v 



(6) 



We restrict our attention to the v > min{^,^ c } cases, for 
which (3 > 1/2, as observed by XRT for the Swift X-ray 
afterglows. From equations (3)-(5), it is easy to obtain the 
synchrotron light-curve decay F„ oc t~ a : 



(Sa) a(vi <v< fc)-|8 = 



(Sc) a(v > v c ) 



(ja) a{vi <v< ^c)~-/3 



6-s 
8- 2s 



(jc) a (v>v c )-^=f-± 




-1/2 


S — > 





s = 


1/2 


s = 


3/2 


s — 


1/2 s 




3/4 s 


= 


1 s 


= 2 


3/2 s 


= 3 


S — ► 


— oo 



(7) 



(8) 



electrons radiating synchrotron emission at the observing 
frequency are losing energy adiabatically) , while the models 
with v c < v will be called Sc , jc , and Jc (where the letter 
"c" shows that the electrons radiating at v arc cooling radia- 
tively). Note from equations (8), (11), and (12) that for the 
Sc and J models, the index a is independent of the medium 
structure, hence the type of CBM cannot be determined for 
these models. 



2.2 Inverse-Compton Dominated Electron 
Cooling 

In the derivation of equation (3) we have ignored a multi- 
plicative factor (Y + 1)~ 2 (where Y is the Compton parame- 
ter) in the expression of v c . Therefore equation (3) is valid if 
Y < 1 (i.e. the radiative cooling of the electrons emitting at 
v c is synchrotron-dominated) or if Y constant (which corre- 
sponds to the v c < Vi case, where the Y parameter depends 
only on the ratio of the electron and magnetic field ener- 
gies). If Y > 1 and Vi < v c , the decrease of the Compton 
parameter with time leads to a faster increase or a slower 
decrease of v c than given in equation (3) and to a slower 
decay of the afterglow emission at v > v c . This case is most 
likely relevant for the 2 nd X-ray afterglow phase, between 
the flattening and steepening times £f and ts, when the 
X-ray light-curve may exhibit a slower decay than that re- 
sulting from equations (8), (10), and (12). The equations for 
the afterglow light-curve at v > v c for the (Y > 1, Vi < v c ) 
case, derived by Panaitescu & Kumar (2001), lead to: 



(9) 



.(10) 



(Sc,s = 0) a--p 



From the above equations, it can be seen that the pas- 
sage of the cooling frequency through the observing band 
steepens the afterglow decay by Aa = |4 — 3s|/(16 — 4s), 
which is at most 1/4, in addition to softening the spectrum 
by A/3 = 1/2 for s < 4/3 or hardening it by A/3 = -1/2 
for s > 4/3. The representative values chosen for s these 
equations show that the observable quantity a — 1.5/3 has a 
stronger dependence on the CBM structure for s 3 (winds) 
than for s ~ 0. The case s = 3 should be taken only as the 
s — > 3 limit; for s = 3 the outflow deceleration is not a 
power-law in the observer time, instead T oc 1/Vlnt. 

For the J model, the (a, 0) closure relation is: 



(Ja) a(vi < v < Vc) - 2/3 
(Jc) a(v > v c ) - 2/3 = . 



(11) 
(12) 



Equations (8), (10), and (12) are valid whatever is the 
location of the injection frequency. However, there are fur- 
ther constraints for the applicability of the v c < v < Vi case: 
(3=1/2 for all models and a = 1/4, 1 for the Sand J models 
respectively. 

The models S , j , and J with v < v c will be desig- 
nated as Sa , ja , and Ja (the letter "a" indicating that the 



-1/(4-2/3) p<3 
-1 p > 3 



(13) 



Ucs = 0) a -\p={%™*-W HI (15) 

0,-2) { a" P<3 (16) 

For simplicity, the results in equations (13)-(16) are given 
for the two most likely types of CBM structure - s = and 
s = 2 - and not for any s. For the J model, a is quasi- 
independent on the stratification of the CBM: 



( Jc) a - 2/3 = 



(l-/3)/(2-/3) p<3 
-1 p > 3 



(17) 



The results given for p < 3 in equations (13)-(17) are 
valid also for p < 2 as long as the total electron energy is a 
constant fraction of the post-shock energy, which is equiv- 
alent to saying that the high-energy cut-off of the electron 
distribution, which must exist for a finite total electron en- 
ergy, has the same evolution as the minimum electron en- 
ergy, 7i oc r. 

The above equations show that the passage of the cool- 
ing frequency through the observing domain slows the after- 
glow decay by Aa > —1/4 for s = and Aa > —5/4. 
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2.3 Radiative Afterglows 

The temporal evolutions given in equations (3) and (5) were 
derived under the assumption of an adiabatic blast-wave. 
If the electron fractional energy is around 50 percent and 
if the electrons cool radiatively (v c < z/j), then radiative 
losses become important. In this case the afterglow emission 
decays faster than for an adiabatic GRB remnant, given the 
stronger deceleration, therefore radiative blast-waves should 
be of importance for the fast decaying, very early Swift X- 
ray afterglows. 

From i) the dynamics of a fully radiative blast-wave 
(FM = const, where M oc nR 3 oc R 3 ~ s is the mass of 
the swept-up CBM) and using ii) the scalings for the spec- 
tral characteristics (vi, c oc 7 ic BY, where B oc Tn 1/2 is the 
post-shock magnetic field strength, 7^ oc V is the electron 
energy, and 7 C oc T 3 r~ 2 B~ 3 is the energy of the electrons 
whose radiative cooling timescale is equal to the dynamical 
timescale; F p oc YBM for the S model and F p oc Y 3 BM for 
the j model), and Hi) the relation between the observer time 
and blast-wave radius r oc T 2 t, the following evolutions of 
the spectral characteristics can be derived: 

(RS, Rj) - n _ = - II _ 5 -, = (is) 



dint 



(RS) 



dlnFp 
dint 



6 - 



14 -4s 



, D ., dlnF p 18 -5s 

(Rj) 



dint 



14 -4s 



(19) 



(20) 



Hereafter, radiative afterglows will be indicated with the 
letter "R" preceding the specific model. Note from equa- 
tion (18) that, just as for an adiabatic afterglow, the cooling 
frequency increases for s < 4/3 and decreases for s > 4/3. 

The light-curve decay indices resulting from equations 
(6) and (18)-(20) are: 



(RSc) a(y c < v < vi) = 

(RSc) a(v c < Vi < v) = 

(Rjc) a(v c < v < Vi) = 

(Rjc) Oe(v c < Vi < v) = 



16 - 5s 
28 - 8s 

24/3 - 4 - s(7/3 - 1) 

14 -4s 
40 - 13s 



28 -8s 



24/3 + 8 - s(7f3 + 3) 
14 -4s 



(21) 
(22) 
(23) 
(24) 



The condition v c < z/j required by radiative dynam- 
ics guarantees that the Compton parameter Y is constant, 
hence there are no further complications with the inverse 
Compton-dominated electron cooling, as it was the case for 
an adiabatic blast- wave. Given that, in the J model, the jet 
Lorentz factor decreases exponentially with radius (Rhoads 
1999), the dynamics and light-curves of a radiative jet should 
be close to those for an adiabatic jet (eqs.[ll] and [12]). 

2.4 Structured Outflows 

There are two other factors which can alter the afterglow 
decay index a. One is that the relativistic outflow can be 
endowed with an angular structure, where the ejecta kinetic 
energy per solid angle, dE/dfi, is not constant (Meszaros, 
Rees & Wijers 1998). The light-curve decay indices for an 
axially-symmetric outflow with a power-law structure 



(25) 



where the angle 9 is measured from the symmetry axis 
(which, for simplicity, is assumed to be also the direction 
toward the observer) are given in Meszaros, Rees & Wijers 
(1998) and Panaitescu & Kumar (2003). In this work, re- 
course to a structured outflow will be made only to explain 
afterglow decays which are slower than that expected for the 

5 model. Evidently, such structured outflows require q > 0. 
If the slow X-ray decay is preceded by a faster fall-off, then 
the index q changes to q < close to the outflow axis, cor- 
responding to the j or J models. If the slow X-ray decay 
is followed by a steepening, then, going away from the out- 
flow axis, the index q changes to either q — (if the steeper 
decay is accommodated by the S model) or to q < (if 
that steeper decay can be explained with the j and J mod- 
els). Therefore, in the most general case, where the X-ray 
light-curve exhibits a sharp decay followed by a slow fall- 
off and then a steeper dimming, the outflow should have a 
bright spot moving toward the observer, surrounded by a 
dim envelope (so that a steep decay is obtained when the 
spot edge becomes visible to the observer) , which is embed- 
ded in a more energetic outer outflow (yielding the slower 
decay) , whose collimation leads to the late steepening when 
the outflow boundary becomes visible. 

The decay index for the synchrotron emission from a 
structured outflow can be derived as described in Panaitescu 

6 Kumar (2003). For a power- law radial structure of the 
CBM and angular structure of the outflow, we obtain: 

^< V <u e) -^=l^f^±^ (26) 
«(, > ,c) ^ \P ^ _ 4- a + 0.W- + ^ + 4- 8 ] (27) 
The above results are valid for 



q > q ■ 



-2(4- 



[3 + 2/3-0.5s(/3 + 3)]- 1 v, < v < v c 
[3 + 2/3-0.5s(/3 + l)n 1 Vi,v c <v 



because for q < q < the emission from the outflow axis 
(6 = 0), where the energy per solid angle would formally 
diverge, becomes dominant and sets another light-curve de- 
cay index. From equations (26) and (27) it follows that, for 
a given CBM structure, the slowest decay that a structured 
outflow can produce is that obtained in the q — > 00 limit: 



-3 + 0.5(/3 + 3)s v t <v < v c 
-2 + 0.5(/3 + l)s Vi,v c <v 



(29) 



Hence, for a homogeneous medium (s = 0), the light-curve 
from a structured outflow could rise (a m i„ < 0). Evidently, 
the structured outflow model can be at work only if the 
above decay is slower than that observed, the condition a > 
otmin leading to a constraint on the CBM structure: 



s < s„ 




(30) 



Equations (26) and (27) give the outflow structural pa- 
rameter which accommodates the observed light-curve index 
a and spectral slope (3: 
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(4 - s)(6/3 - 4a + 3) + 5s - 12 

2a + 6- ((3 + 3)s 
(4-s)(6/3-4q) + 2s-8 
2a + 4- {/3 + l)s 



Vi, v c <v 



.(31) 



2.5 Energy Injection 

Another process which can reduce the afterglow dimming 
rate is the injection of energy in the blast-wave (Paczyhski 
1998, Rees & Meszaros 1998) by means of some ejecta which 
were ejected later than the GRB ejecta (a long-lived engine) 
or at the same time but with a smaller Lorentz factor, thus 
reaching the decelerating GRB ejecta during the afterglow 
phase (a short-lived engine). A delayed injection of energy 
into the afterglow can be due to the absorption of the dipole 
electromagnetic radiation emitted by a millisecond pulsar 
(Dai & Lu 1998, Zhang & Meszaros 2001) if such a pulsar 
was formed. 

The addition of energy in the blast-wave mitigates its 
deceleration and, implicitly, the afterglow decay rate. Rees 
& Meszaros (1998) have derived the decay index a for an 
energy injection that is a power-law in the ejecta Lorentz 
factor. The expressions for the index a for an energy injec- 
tion which is a power-law in the observer time, 



Ei(< t) tx f 



(32) 



are given in eqs.(23), (24), and (30) of Panaitescu & Kumar 
(2004). From those equations, it follows that energy injec- 
tion reduces the light-curve decay indices given in equations 
(2.1)-(12) by 




(j-EI) Aa = e- i f 



(J -EI) Aa=-e- 



jy(/3 + 2) + g 2 _| g Vi <v<v 
+ + i>i,Vc<v 



(3 + 2 Vi < v < v c 
(3 + 1 Vi,v c <v 



(34) 



(35) 



for the adiabatic S , j , and J models. 

Lastly, all the decay indices given in the above equations 
were derived assuming that the microphysical parameters 
which determine the spectral characteristics (z/», i/ c , F p ) and 
the continuum slope (/3), i.e. the parameters for the typical 

post-shock electron energy & magnetic field strength^ and 
the power-law index p of the electron distribution with en- 
ergy, are constant. This possibility is not investigated in this 
work. 



3 MODELS FOR SWIFT X-RAY 
AFTERGLOWS 

As described in the Introduction, the Swift X-ray afterglows 
exhibit three phases: the l si phase, lasting until tF ~ 300 
s, is characterized by a sharp decay, the 2 nd phase, lasting 



t Yost et al. (2003) have shown that a decrease of the parameter 
for the magnetic field energy slower than i^ 1 / 2 or an increase 
slower than t 3 / 4 are allowed for several afterglows 



until ts ~ 10 3 ' 5 — 10 5 s, is marked by a much slower fall- 
off, while in the 3 rd phase, the X-ray light-curve displays a 
faster decay. The light-curve decay indices a and the spec- 
tral slopes (3 are listed for each phase in table 1. The closure 
relations between a and (3 presented in section 2 provide ei- 
ther a criterion for distinguishing among the various models 
that can accommodate the observed afterglow properties or 
allow the determination of the CBM structure. Since s < 3 
is required for a decelerating blast- wave, this also serves as 
a test of the various models. 

Table 1 lists the models for which the closure relations 
given in section 2 between the light-curve decay index a and 
spectral slope (3 are satisfied within la, for each afterglow 
decay phase. To find a model for the entire afterglow, these 
piece-wise models must now be put together in a sequence 
that makes sense and is not contrived. The criteria by which 
we construct a model for the entire X-ray afterglow are: 

i) models relying on coincidences to accommodate two adja- 
cent X-ray phases are excluded, i.e. only one factor (cooling 
frequency passage, change of CBM structure, region of non- 
monotonic variation in the energy per solid angle becoming 
visible, beginning/cessation of energy injection) at a time is 
employed to explain a variation of the X-ray decay index, 

ii) radiative outflows can evolve into adiabatic ones, but not 
the other way around, 

iii) any of the three dynamical models (S , j , J ) can be 
followed by the same model, but only the S model can be 
followed by the j and J models, allowing for a collimated 
outflow, spreading or non-spreading, whose edge becomes 
visible to the observer, 

iv) the evolution of the cooling frequency v c required to join 
two models at If or ts must be compatible with the CBM 
structural index s, i.e. v c can increase only if s < 4/3 and 
can decrease only if s > 4/3 (modulo the effect of a decreas- 
ing Compton parameter when electron cooling is dominated 
by inverse Compton scatterings) . 

A structured outflow or energy injection are invoked 
only when the X-ray decay during the 2 nd phase (tp < t < 
ts) is too slow to be explained by the S , j , and J models, 
for the CBM structured required to accommodate the X-ray 
decay preceding (t < t F ) or following (t > ts) this phase: 

v) for the structured outflow model, a working condition is 
that the slowest decay that it would yield (eq.[29]), given the 
CBM structure which explains the X-ray emission at t < t F 
or at t > ts, is slower than that observed. A structured out- 
flow for the 2 nd afterglow phase cannot be preceded by the 
J model, as the existence of an outflow outside the jet would 
prevent its lateral spreading, 

vi) for the energy injection model, we determine from equa- 
tions (33)-(35) the index e (eq. [32]) which reconciles the 
slow X-ray decay with the spectral slope, for the model (S , 
j , or J ) and CBM structure which accommodates the X- 
ray emission before (t < tp) or after (t > ts) the energy 
injection episode. The ratio Ei/Eo of the total injected en- 
ergy Ei to the energy Eo existing in the blast-wave prior 
to the energy injection episode is (t ff/t on ) e where t on is 
the light-curve flattening time tF or the epoch of the first 
measurement (if no flattening was observed) and t a f / is the 
light-curve steepening time ts or the epoch of the last mea- 
surement (if no steepening was observed). Then a test of the 
energy injection model can be done if it is assumed that the 
pre-injection energy Eo is comparable to the GRB 15-350 
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Table 1. Models that accommodate the light-curve power-law decay index (F x oc t a ) and continuum power-law slope (F v tx v I 3 ) 
measured by Swift for the early (~ 10 2 s), mid (10 3 — 10 4 s), and late (~ 10 5 s) X-ray afterglow emission 





l si phase: t < tp 


- Steep Decay 


2 nd phase: t F <t <t s - Slow Decay 


3 rd phase: t s < t 


- Fast Decay 


GRB 




Pi 


Model 






Model 






Model 


050126" 


2.7±.2 


1.3±.2 


Sa,ja,Rjc,Jc 


0.6±.l 


~ 1.3 


none 








050128 6 








0.3±.l 


0.6±.l 


Sc (jc,RSc,Rjc) 


1.3±.2 


0.8±.l 


Sa,jc,Jc (RSc,Rjc) 


050219A C 


3.2±.2 


1.1±.2 


Saja, Ja 


0.8±.l 


~ 1.1 


Sc (RSc,Rjc) 








050315° 


3.3±.2 


1.3±.2 


Saja 


0.7±.l 


0.9±.2 


Scjc (RSc.Rjc) 


1.6-2.6 


1.0±.l 


Sajajc,Ja,Jc (RSc,Rjc) 


050318 a 








l.Oi.l 


1.0±.2 


Sa,Scjc,RSc (Rjc) 


1.4±.3 


1.2±.3 


Sa,jajc,Jc (RSc,Rjc) 


050319" 


3.0±.2 d 


1.9±.2 


Sajc, RSc, Rjc 


0.5±.l 


0.8±.l 


Scjc, RSc, Rjc 


1.2±.l 


0.7±.2 


Sajc,Jc (Rjc) 


050401" 








0.5±.l 


l.li.l 


none 


1.6±.l 


l.li.l 


Sajc (RSc,Rjc) 


050408" 








0.7±.l 


1.1±.2 


none 


1.5±.2 


~ 1.1 


Sajc (RSc, Rjc) 


050505" 








0.7±.2 


l.Oi.l 


Scjc (RSc.Rjc) 


2.5±.4 


0.9±.l 


Saja,Ja (Rjc) 



Refs. for a and f3 - " Chincarini et al. (2005), b Campana ct al. (2005), c Tagliaferri ct al. (2005) 
d for t = at the beginning of the second GRB peak (Barthelmy et al. 2005) 

Model coding - S: spherical outflow, j: non-spreading jet, J: sideways spreading jet, R: radiative afterglow 

(for models given in parentheses the outflow is less likely to be radiative at ~ 1 day after the burst, or require a 
wind with a radial profile close to r , for which the analytical results given in section 2 are only approximative) 
a: v x < v c (X-ray emitting electrons are cooling adiabatically) 
c: v c < v x (X-ray emitting electrons are cooling radiativcly) 



keV output. For this test, we calculate the injected energy 
Ei and require for the S model that the outflow kinetic 
energy contained within 10° (which is sufficiently wide to 
resemble a spherical outflow until 1 day after the burst) 
does not exceed 10 53 ergs. However, this test is based on the 
assumption that the ejecta producing the burst and the be- 
ginning of the X-ray afterglow emission are the same. This 
does not have to be the case, as it is possible that the GRB 
emission arises from internals shocks in the entire outflow 
and yet only its leading edge drives the forward shock and 
radiates at the beginning of the X-ray afterglow. 

We allow a variable index s as we do not know what 
are the properties of the winds expelled by massive stars 
in the last 1,000 years before they explode, hence we do 
not know what is the density structure of the CBM within 
the first parsec, where the afterglow emission is produced. 
Variations in the GRB progenitor's mass-loss rate and wind 
speed could lead to a CBM with a structure different than 
the r~ 2 profile expected for a constant speed, constant mass- 
loss rate, and to interactions between winds that could form 
shells of higher density. 

By applying the above criteria, we arrive at the models 
given in Table 2. Below we discuss in some detail the 9 Swift 
X-ray afterglows and the models that accommodates them. 

050126. The XRT light-curve of this afterglow exhibits 
a steep fall-off until tF <) 300 s, followed by a slower decay. 
Tagliaferri et al. (2005) have shown that extrapolation of the 
15-350 keV BAT emission to the 0.2-10 keV XRT band is 
dimmer at 100 s than the observed XRT flux. Furthermore, 
the XRT spectrum = 1.26 ± 0.22) during the I s * phase 
is softer than the BAT spectrum (/3 7 = 0.32 ± 0.18), hence 
the early X-ray afterglow is not the large-angle GRB emis- 
sion and must be attributed to the forward shock. If there 
is no spectral evolution (P2 = Pi) across tF, as indicated by 
Tagliaferri et al. (2005), then the slow X-ray decay of the 2 nd 
phase cannot be explained by a change in the structure of the 
CBM medium for any of the models (Sa , ja , Rjc , Jc ) which 
accommodate the 1 st phase. Conversely, if the CBM struc- 
ture does not change across tF, then the slower decay at the 



2 phase requires a substantial hardening of the spectrum, 
corresponding to a rising one (P2 < 0) for the models Sa , 
ja , and Rjc , or one with (3? = 0.20 ± 0.25 for the Jc model, 
both of which are inconsistent with the XRT observations. 
Furthermore, for the possible models for the I s afterglow 
phase, the passage of the cooling frequency through the X- 
ray band can only steepen the afterglow decay. Hence, the 
most plausible models that can explain the flattening X-ray 
light-curve of 050126 require energy injection or a struc- 
tured outflow. The Sa model with either energy injection or 
a structured outflow does not satisfy conditions v) and vi) 
above. 

050128. Although XRT observations started at 100 s 
after the burst, a steep early decay has not been observed 
(Campana et al. 2005). Its decay steepens at ts ii 10 3 s, 
without a spectral evolution. Of the many possible com- 
binations of models for the 2 nd and 3 rd phases, the most 
plausible is that of a collimated outflow (jc and Jc models) , 
leading to a steepening of the X-ray decay when the bound- 
ary of the jet becomes visible. Another possibility is that 
of non-spreading jet (jc model) which transits from a r~ 3 
wind into a region of increasing density at t,s- We note that 
all these models require a rather hard electron distribution, 
with p J5 1.3. 

050219A. The features of this afterglow are similar to 
those of 050126. It exhibits a fast fall-off until t F ~ 300 s, 
followed by a slower decay. The extrapolation of the 15-350 
keV BAT emission to the 0.2-10 keV XRT band underpre- 
dicts the observed flux at 100 s (Tagliaferri et al. 2005) and 
the X-ray spectral slope (/3i = 1.1±0.2) is much softer than 
that of the burst (/3 7 = — 0.75±0.30), hence the rapid, early 
fall-off of the 050219A X-ray afterglow is not the GRB large- 
angle emission. Just as for the afterglow 050126, a change 
in the CBM structure cannot explain the X-ray light-curve 
flattening. If the CBM structure is considered unchanged 
across tF, then the slowing of the X-ray decay would re- 
quire a rising spectrum (/?2 < 0) for the 2 nd phase, which is 
inconsistent with the XRT observations. Because all models 
for the I s afterglow phase require that the cooling frequency 
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Table 2. Possible models for the afterglow X-ray phases of Table 1 



t < tp - Steep Decay tp < t < tg - Slow Decay tg < t - Fast Decay 

GRB Model s Model s e q Model s p 







(1) 




(1) 


(2) 


(3) 




(1) 


(4) 


050126 ja 




< 2.5 


EI 




1.0-1.4 








3.5±.4 


ja 




< 2.5 


SO 


< 1.7 




> 0.7 






3.5±.4 


Rjc 




< 3 


EI 




1.0-1.9 








2.5±.4 


Kjc 




< O 




< Z.O 




> 1.4 






Z.Ozt.4 


Jc 




< 3 


EI 




1.4 








2.5±.4 








OC 


< O 






jc(Jc) 


< l.l\6 ) 


1 o_l_ 1 








jc 


3 






jc 


< 


1.2±.l 


050219A ja 




2.4-3 


EI 




1.7-2.3 








3.2±.4 


Ja 




< 3 


EI 




1.2 








3.2±.4 


050315 LA- 






Sc 


< 3 






jc(Jc) 


< 2.1(3) 


1.9±.l 


050318 






Sc 


< 3 






jc(Jc) 


< 1.9(3) 


2.1±.l 








c 

Sa 


< —3 






ba 


2 


3.0±.3 








jc 


2.6-3 






jc 


— ► —oo 


2.1i.l 


050319 LA- 


GRB 




Sc 


< 3 






jc(Jc) 


< 2.1(3) 


1.5±.2 


LA- 


-GRB 




jc 


3 






jc 





1.5±.2 


050401 






EI 




0.7 




Sa 


-1-0.6 


3.1±.l 








SO 


< 1.7 




1.7-3.5 


Sa 


-1-0.6 


3.1±.l 








SO 


< 2.4 




> 1.1 


RSc 


< 3 


2.1±.l 








EI 




0.7 




jc 


1.6-2.4 


2.1±.l 








SO 


< 2.4 




> 3.2 


jc 


1.6-2.4 


2.1±.l 








EI 




0.6 




Rjc 


2.7-2.9 


2.1±.l 


050408 






EI 




0.4-0.6 




Sa 


< 1.1 


3.3±.4 








SO 


< 1.8 




> 0.5 


Sa 


< 1.1 


3.3±.4 








EI 




0.4-0.6 




jc 


0.7-3 


2.3±.4 








SO 


< 2.5 




> 0.8 


jc 


0.7-3 


2.3±.4 


050505 






EI 




> 0.9 




ja/Ja 


< 3 


2.8±.2 








SO 


< 1.9 




> 0.9 


ja 


-2—3 


2.8±.2 



Model coding - EI: energy injection; SO: structured outflow; LA-GRB: large-angle (9 > 1/r) GRB emission 

(1) : exponent of radial density profile (eq.[2]); for the SO model, the upper limit on s is that resulting from equation (30) 

(2) : exponent of the energy injection law (eq.[32]) obtained from eqs.(33)-(35) for the index s required at t < tp or at t > tg 

(3) : exponent of the angular distribution of the energy per solid angle (eq.[25]) obtained from equation (31) for the index s required 
at t < tp or at t > tg 

(4) : exponent of the power-law distribution of electrons with energy (eq.[l]); this value is for all X-ray phases except LA-GRB 



is above the X-ray domain, its passage is either impossible 
or it would steepen the light-curve decay. Consequently, the 
slower decay observed for the X-ray afterglow 050219A after 
tp requires either energy injection or a structured outflow. 
Condition v) is not satisfied by either the Sa and ja models 
and a structured outflow, while the Sa model with energy 
injection requires too much energy. 

050315. The X-ray emission exhibits a flattening at 
tp ^ 10 3 s, accompanied by a hardening of the spectrum 
(P'2 — f3i = —0.41 ± 0.21), and followed by a steeper decay 
after tg i; 10 s, across which there is no spectral evolution. 
Barthelmy et al. (2005) have shown that the early, steep 
fall-off is consistent with the large-angle GRB emission: the 
extrapolation of the 15-350 keV BAT emission to the 0.2-10 
keV XRT band matches the XRT flux measured at 100 s, the 
X-ray spectral slope (f3i = 1.34±0.15) is comparable to that 
of the burst (/3 7 = 1.18 ± 0.11), and the X-ray decay index 
(qi = 3.35 ± 0.19) is close to the expected value (2 + /3 7 = 
3.18±0.11). The steepening at tg can be easily understood as 



due to a collimated outflow (the j or J models). A radiative 
non-spreading jet interacting with s < 3 CBM could also 
accommodate the steepening, if the CBM is a wind, but it 
is less likely that the radiative phase could last until later 
than 1 day after the burst. 

050318. Because XRT observations started at 1 h 
after the burst, the fast decay phase may have been missed. 
A steepening of the X-ray light-curve decay occurs at tg ~ 
3 x 10 4 s without a spectral evolution. This steepening can 
be due to seeing the boundary of a jet (spreading or not). 
There are other possible models that can accommodate the 
steepening, all involving a variation in the CBM structural 
index s. They are the Sa outflow exiting a shell of a sharply 
increasing density and entering a r~ 2 wind and jc outflow 
transiting from a r~ 3 wind to a shell with sharply increasing 
density at is. 

050319. This afterglow is similar to 050315, the hard- 
ening of the X-ray spectrum across the light-curve flatten- 
ing, which occurs at tp ~ 400s, being stronger. Barthelmy 
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et al. (2005) have shown that the BAT GRB emission ex- 
trapolated to the XRT band matches the X-ray flux mea- 
sured at 200 s. If the origin of time for the X-ray emis- 
sion is set at the beginning of the second (and last) GRB 
pulse, then the decay index (ai = 3.0 ± 0.2) of the early 
X-ray emission is consistent with the expectations for the 
large-angle GRB emission (2 + /3 7 = 3.13 ± 0.28). However, 
the early X-ray spectrum (/3i = 1.94 ± 0.20) is rather soft 
compared to that of the burst. On the other hand, the sub- 
stantial hardening of the X-ray spectrum across tF, with 
02 — f3i = —1.15 ± 0.23, exceeds that which the passage of 
the cooling frequency through the observing band can pro- 
duce (P2 — Pi — —0.5), suggesting that the X-ray emissions 
during the 1 st and 2 nd afterglow phases arise from different 
mechanisms. We note that the X-ray light-curve for both 
phases may be explained in the structured outflow frame- 
work if we make the ad-hoc assumption that the spectrum 
of the spot emission (dominating the afterglow flux before 
ti?) is softer than that from the surrounding outflow (which 
overtakes the spot emission after tF)- 

The steepening of the X-ray light-curve at ts ~ 3 x 10 4 
s can be explained by seeing the edge of a jet (spreading or 
not), or with the jc model and a CBM structure changing 
from a r~ 3 wind to a homogeneous medium at ts- All these 
models require a hard electron distribution, with p < 1.7. 

050401. Although the XRT observations started 100 
s after the burst, a steeply falling-off phase was not seen. 
Until ts — 5000 s, it exhibits a decay so slow that it can- 
not be explained without energy injection or a structured 
outflow. The RSc model with energy injection requires too 
much energy, while the structured outflow does not satisfy 
condition u) for the Rjc model. Then the steepening of the 
X-ray light-curve at ts can be understood either as resulting 
from the cessation of the energy injection or from seeing the 
outflow boundary. In the latter case, the light-curve decay 
should be steeper than for the S model and slower than for 
the j model. That the steeper decay after ts can be accom- 
modated by either the S and j models (Table 1) supports 
a structured outflow as the source of the X-ray light-curve 
steepening. 

050408. This afterglow is very similar to 050401, except 
that the steepening occurs later, at ts ~ 10 5 s. Its light- 
curve decay before ts is also too slow and requires an energy 
injection episode or a structured outflow. The X-ray spectral 
slope after ts is not known, but if we assume that there is 
no spectral evolution across ts (as is the case for all other 
afterglows), then the light-curve decay index and spectral 
slope measured after ts can be accommodated by the Sa and 
jc models. The RSc and Rjc models are also allowed, though 
it is unlikely that the radiative phase lasts until days after 
the burst. 

050505. This afterglow is similar to 050318, its X-ray 
light-curve steepening at ts = 4 x 10 4 s without a spec- 
tral evolution. However, in contrast to 050318, the spectral 
slopes and decay indices before and after ts cannot be rec- 
onciled within any model other than Rjc , even if we al- 
low for a varying CBM structure. Besides that the radia- 
tive phase is unlikely to last until 1 day after the burst, the 
Rjc model requires a r~ 3 CBM profile, for which the clo- 
sure equations given in section 2.3 are not accurate. Hence, 
it seems more plausible that the slow decay of this after- 
glow before ts is due to energy injection or a structured 



outflow. The Sa model fails to satisfy conditions v) and vi) 
for these two case. As for the afterglows 050401 and 050408, 
the steepening of the X-ray light-curve could then be at- 
tributed to the end of the energy injection or to the outflow 
axis becoming visible to the observer. 



4 DISCUSSION 

As shown in Table 1, the three decay phases of the Swift 
X-ray afterglows can be understood in the following way: 

i) the hardening of the 0.2-10 keV spectrum of the X-ray 
afterglows 050315 and 050319 from t < 400 s (when a fast 
decaying X-ray emission is observed) to t > 400 s (when the 
X-ray light-curve exhibits a slow decay) indicates that the 
early, fast-falling off X-ray emission arises from a different 
mechanism than the rest of the afterglow. Barthelmy et al. 
(2005) have argued that this mechanism is the same as for 
the GRB emission. The results of Tagliaferri et al. (2005) 
suggest that this explanation does not work well for the 
afterglows 050126 and 05018. In these cases, the steep X- 
ray decay require a very narrow outflow whose edge is seen 
as early as 100 s after the burst, the following, slower decay 
phase being explained by energy injection in the blast-wave 
or by a rather contrived (see below) angular structure of the 
blast- wave, 

ii) the X-ray decay measured until Hi, 0.3d, and 0.5d for the 
afterglows 050401, 050408, and 050505, respectively, is too 
slow to be explained by the simplest blast- wave model. Such 
a slow decay can be produced by an outflow endowed with 
angular structure or by a continuous injection of energy in 
the forward shock, 

in) the steepening of the X-ray decay observed at lh-2d for 
the afterglows 050128, 050315, 050318, and 050319, i.e. at 
a time comparable to that of the steepening of the optical 
decay of many pre- Swift afterglows, can be explained by 
seeing the edge of a jet. For the remaining 5 afterglows, 
whose pre-break X-ray decay requires energy injection, the 
steepening can be attributed to the cessation of injection. 

In the large-angle GRB emission model for the early, 
fast falling-off phase, the X-ray emission arises from the 
same mechanism as the GRB itself, but arrives at observer 
later because it comes from the shocked gas moving slightly 
off the direction toward the observer. For this model to be at 
work, three conditions must be satisfied. First, the 15-350 
keV GRB emission extrapolated to the 0.2-10 keV X-ray 
band, under the assumption that the power-law burst spec- 
trum F v oc v~^~< extends unbroken down to 0.2 keV, should 
match or exceed the X-ray flux at the first epoch of ob- 
servations. Second, the spectral slope of the early afterglow 
should be the same as that of the GRB. Third, the X-ray 
light-curve decay index should be equal to 2 + (3 1 (Kumar & 
Panaitescu 2000). For completeness, we present here a short 
derivation of this result. If the GR emission stops suddenly 
at some radius r and blast-wave Lorentz factor T, then the 
received flux is F„(t) oc J^, (d£/dt)£> 2 where T' u , oc v'~ 9 is 
the outflow comoving frame surface-brightness at frequency 
v' — v/T>, dE = 2-kt 2 8A9 is the elementary area whose ra- 
diation is received over an observer time dt, 6 is the an- 
gle (measured from the direction toward the observer) of 
the fluid element from which radiation is received at time 
t = r9 2 / 2 (hence dt oc 6 AO), V = 2 /(TO 2 ) oc i" 1 is the 
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relativistic Doppler factor, and the expressions for t(9) and 
V(9) have been derived for 8 3> r _1 , i.e. for the large-angle 
emission. The last factor T> 2 in the expression of F v accounts 
for the beaming of radiation from a relativistic source. After 
substitutions, one obtains that F v oc i/~ f3 t~ 2 ~ f3 . 

Barthclmy et al. (2005) found that the above three 
conditions for the large-angle GRB emission as the source 
of the very early, fast X-ray decay are met for the after- 
glows 050315 and 050319. For two other afterglows, 050421 
(Sakamoto et al. 2005, Godet et al. 2005) and 050713B (Par- 
sons et al. 2005, Page et al. 2005), we find that their fast X- 
ray decays cannot be reconciled with their hard X-ray con- 
tinua by any of the blast-wave models considered here, but 
they do satisfy the last two conditions above for the large- 
angle GRB emission interpretation. Tagliaferri et al. (2005) 
showed that the early X-ray emissions of the afterglows 
050126 and 050219A are brighter than the GRB extrapo- 
lated fluxes and softer than the burst emission, therefore 
their early X-ray afterglows cannot be immediately identi- 
fied with the large-angle GRB emission. Kumar et al. (2005) 
discuss the conditions under which the fast X-ray decay of 
these last two afterglows can be reconciled with the large- 
angle GRB emission. 

The outflow structure required to explain a light-curve 
flattening followed by a steepening must contain a bright 
spot (moving toward the observer) surrounded by a dimmer 
envelope where the ejecta have a lower energy per solid an- 
gle dE/dfi, so that the emission from the spot exhibits a 
fast decay after its boundary becomes visible. Further, the 
envelope should be embedded in a wider outflow whose emis- 
sion overtakes that of the spot when the blast-wave Lorentz 
factor has decreased sufficiently. To explain the light-curve 
flattening, the dE/dfi in this wider outflow should rise away 
from the spot as 6 |1/ ' 2 to 9 3 , where 6 is the angle measured 
from the outflow's symmetry axis. Finally, to explain the 
light-curve steepening, the dE/dfi should stop increasing 
with angle (for the S model) or peak and then decrease (for 
the j and J models). The decrease could be gradual, with the 
X-ray light-curve steepening occurring when the fluid at the 
peak of dE/dfi becomes visible and the outer outflow con- 
tributing to the post-break emission (this is the light-curve 
break mechanism proposed by Rossi, Lazzati & Rees 2002). 
If the decrease of dE/dfi is sharp, then the post-break X-ray 
light-curve decay will be faster, particularly if the outflow 
undergoes lateral spreading (this is the light-curve break 
mechanism proposed by Rhoads 1999). 

In the energy injection model, the forward shock energy 
increases due to some relativistic ejecta which catch up with 
the decelerating blast-wave (Rees & Meszaros 1998). The en- 
ergy injection reduces the blast-wave deceleration rate and 
mitigates the decay of the afterglow emission. In principle, 
during the slow X-ray decay phase, there could be an en- 
ergy injection for all afterglows considered here; Table 2 lists 
only the cases when it is required. As shown in Table 2, we 
find that, to explain the slow phase of the X-ray afterglow 
decay, the blast-wave energy should increase with observer 
time faster than t B and slower than t 1 ' 5 . The arrival of 
new ejecta at the forward shock could be due either to the 
spread in the Lorentz factor of ejecta released simultaneously 
(short-lived engine) or to a long-lived GRB engine, releasing 
a relativistic outflow for a source-frame duration compara- 
ble to the observer-frame duration of the slow X-ray decay 



phase. In the former case, the ejecta-forward-shock contrast 
Lorentz factor is [(4 — s)/(l + e)] 1 ^ 2 ^ 2 (Panaitescu & Ku- 
mar 2004), where e is defined by equation (32), while in the 
latter case the Lorentz factor ratio can be much larger. Con- 
sequently, we expect that, for a short-lived engine, the re- 
verse shock propagating in the incoming ejecta is only mildly 
relativistic and radiates mostly at radio frequencies while for 
a long-lived engine the reverse shock could be very relativis- 
tic and radiate in the infrared-optical. 

As mentioned above, the early, fast decay of the X-ray 
afterglows 050126 and 050219A cannot be readily identified 
with the large-angle GRB emission and could be attributed 
to the forward shock. The fast X-ray decay displayed by 
these afterglows at ~ 400 s can then be explained only within 
the structured outflow and energy injection models. In ad- 
dition, two other afterglows, not considered in this work, 
050712 (Grupe et al. 2005) and 050713A (Morris et al. 2005), 
exhibit an X-ray decay which is too slow and incompatible 
with the reported X-ray spectral slope, both requiring either 
an energy injection or a structured outflow. 

For the afterglows 050128, 050318, and 050319, we find 
that a change in the circumburst density profile provides 
an alternate model to structured outflows and jets for the 
steepening of the X-ray decay observed by XRT after 0.1 d. 
For all three afterglows, the changing external density corre- 
sponds to a transition from a r -3 wind to a region of uniform 
or increasing density. The r -3 density structure requires a 
time-varying mass-loss rate and/or speed of the wind of the 
massive star GRB progenitor, while the uniform or increas- 
ing density shell could result from the internal interactions 
in a variable wind (e.g. Ramirez- Ruiz et al. 2005). The self- 
similar solutions of Chevalier & Imamura (1983) for wind- 
wind interactions indicate that a uniform shell results from 
a substantial decrease of the star's mass-loss rate accom- 
panied by a large increase in the wind speed. These major 
changes in the wind properties would have to occur ~ 1, 000 
yrs before the GRB explosion, if the radius where the r -3 
circumburst density profile terminates is that of the forward 
shock at 0.1-1 d. 

To answer the question of how can we distinguish be- 
tween the above three models (energy injection, structured 
outflow, non-monotonic circumburst density) for flattenings 
and steepenings of the afterglow light-curve, we note that, 
if the cooling frequency located between the optical and X- 
ray domains, each of those models yields a specific difference 
Aa c — Aa between the changes Aa c and Aa of the X-ray 
and optical light-curve decay indices. Therefore, to discrimi- 
nate among the few possible models discussed here, it is very 
important to monitor the optical afterglow emission over a 
wide range of times, from minutes to days after the burst. 



5 CONCLUSIONS 

The most often encountered feature resulting from the anal- 
ysis of the nine X-ray afterglows analyzed in this work is the 
existence of a substantial energy injection in the blast-wave 
at hours to 1 d after the burst. This energy injection is nec- 
essary to reconcile the slowness of the X-ray decay with the 
hardness of the X-ray continuum for five out of nine after- 
glows and is possible for the other four as well. The injection 
should increase the blast-wave energy as i 10±0 ' 5 ) t being the 
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observer time, leading to a shock energy which is eventually 
larger by a factor 10-1,000 than that at the beginning of the 
slow X-ray decay phase. 

The exponent e of power-law energy injection identi- 
fied in Table 2 is generally inconsistent with that expected 
for the absorption of the dipole radiation from a millisecond 
pulsar (see also Zhang et al. 2005). In this case, an important 
energy input in the blast-wave is obtained only for the first 
few thousand seconds, when the pulsar electromagnetic lu- 
minosity is constant (Zhang & Meszaros 2001), which leads 
to e = 1. Hence the energy injection must be mediated by 
the arrival of new ejecta at the forward-shock. 

If the GRB engine were so long-lived that only a small 
fraction of the total outflow energy yielded the burst emis- 
sion, then such a large energy injection would imply an even 
higher GRB efficiency than previously inferred (above 30 
percent - Lloyd-Ronning & Zhang 2004) from afterglow en- 
ergetics and would pose a serious issue for the GRB mecha- 
nism (Nousek et al. 2005). However, we cannot yet tell if the 
energy injection lasts for hours because the central engine 
is long-lived or because there is a sufficiently wide distribu- 
tion of the initial Lorentz factor of the ejecta expelled by a 
short-lived engine. In the latter case, the entire outflow could 
be emitting both the GRB and afterglow emission and the 
GRB efficiency remains unchanged. 

If the fast decaying X-ray emission preceding the slow 
X-ray fall-off arises from the forward shock as well, then 
energy injection must start at the end of the steep decay 
phase, i.e. it must be a well defined episode. However, in 
this scenario the outflow must be very tightly collimated to 
yield a fast decay at only 100 seconds after the burst. The 
blast-wave Lorentz factor is T(t = 100s) ~ 100(£ 5 3/no) 1/8 
(E53 being the shock energy in 10 53 ergs and no the circum- 
burst medium particle density in cm" 3 ), thus the jet open- 
ing must be less than = 0.5°. Alternatively, the early 
fast falling-off X-ray emission could be the large-angle emis- 
sion from the burst phase, overshining the forward shock 
emission. This interpretation is favoured by Barthelmy et 
al. (2005) and Hill et al. (2005) for the afterglows 050117, 
050315, 0509319, although for other afterglows (e.g. 050126 
and 050219A - Tagliaferri et al. 2005) the spectral properties 
of the burst and X-ray afterglow emissions are not readily 
consistent with each other. 

For a structured outflow to explain all the three decay 
phases of the Swift X-ray afterglows, the distribution of the 
ejecta kinetic energy with angle must be non-monotonic. 
This is somewhat contrived and inconsistent with the re- 
sults obtained by MacFadyen, Woosley & Heger (2001) from 
simulations of jet propagation in the collapsar model. Thus 
structured outflows do not appear to provide a natural 
explanation for the features of the X-ray afterglow light- 
curves. 

We note that the inferred indices of the power-law elec- 
tron distribution with energy, which are given in Table 
2, range from 1.3 to 2.8. One would have to ignore 4 of 
these 9 Swift afterglows to obtain a unique electron index, 
p — 2.1 ±0.1. This is a puzzling feature of relativistic shocks 
in GRB afterglows: the shock-accelerated electrons do not 
have a universal distribution with energy, a fact which is also 
proven by the wide spread of the high-energy spectral slopes 
of BATSE bursts (A/3 ~ 2.0 in fig. 9 of Preece et al. 2000), 
which is equal to p/2 or (p— l)/2, and the by wide range of 



the optical post-break decay indices of the BeppoSAX after- 
glows (Aa ~ 1.6 in fig. 3 of Stanek et al. 2001 and in fig. 2 
of Zeh, Klose & Kann 2005), which is equal to p. However, 
from the X-ray spectral slope of 15 BeppoSAX afterglows, 
De Pasquale et al. (2005) conclude that the electron index 
p has an universal value of p — 2.4 ± 0.2 (see their fig. 3). 

Acknowledgments. We are indebted to Pawan Kumar for 
his comments on this work. 



REFERENCES 

Akerlof C. ct al. , 1999, Nature, 398, 400 

Barthelmy S. et al. , 2005, ApJL, in press (astro-ph/0511576) 
Blustin A. et al. , 2005, ApJ, accepted (astro-ph/0507515) 
Campana S. etal, 2005, ApJ, 625, L23 
Chevalier L., Imamura J., 1983, ApJ, 270, 554 
Chevalier R., Li Z., 2000, ApJ, 536, 195 

Chincarini G. et al. , 2005, ApJ, submitted (astro- ph/0506453) 
Dai Z., Lu T., 1998, A&A, 333, L87 

De Pasquale M. et al. , 2005, A&A, submitted (astro-ph/0507708) 
Fan Y., Dai Z., Huang Y., Lu T., 2002, Chinese J. A&A, 2, 449 
Fox D. et al. , 2003, ApJ, 586, L5 

Frail D., Waxman E., Kulkarni S., 2000, ApJ, 537, 191 

Frail D. et al. , 2004, ApJ, 600, 828 

Godct O. ct al. , 2005, GCN Circ. 3301 

Grupc D. et al. , 2005, GCN Circ. 3579 

Hill J. ct al. , 2005, ApJ, accepted (astro-ph/0510008) 

in't Zand J. ct al. , 2001, ApJ, 559, 710 

Kulkarni S. et al. , 1999, Nature, 398, 389 

Kumar P., Panaitescu A., 2000, ApJ, 541, L51 

Kumar P., Panaitescu A., 2003, MNRAS, 346, 905 

Kumar P. et al. , 2005, Science, submitted 

Li W., Filippenko A., Chornock R., Jha S., 2003, ApJ, 586, L9 

Lloyd-Ronning N., Zhang B., 2004, ApJ, 613, 477 

Meszaros P., Rees M.J., 1997, ApJ, 476, 232 

Meszaros P., Rees M.J., 1999, MNRAS, 306, L39 

Meszaros P., Rees M.J., Wijers R., 1998, ApJ, 499, 301 

Morris D. ct al. , 2005, GCN Circ. 3606 

Nousek J. et al. , 2005, ApJ, submitted (astro-ph/0508332) 

Paczyriski B., 1998, ApJ, 494, L45 

Page K. et al. , 2005, GCN Circ. 3602 

Parsons A. et al. , 2005, GCN Circ. 3600 

Panaitescu A., Kumar P., 2001, ApJ, 554, 667 

Panaitescu A., Kumar P., 2003, ApJ, 592, 390 

Panaitescu A., Kumar P., 2004, MNRAS, 350, 213 

Pian E. ct al. , 2001, A&A, 372, 456 

Preece R. et al. , 2000, ApJ, 126, S19 

Ramirez-Ruiz E., Garcia-Segura G., Salmonson J., Perez-Rcndon 

B., 2005, ApJ, 631, 435 
Rees M.J., Meszaros P., 1998, ApJ, 496, LI 
Rhoads J., 1999, ApJ, 525, 737 

Rossi E., Lazzati D., Rees M.J., 2002, MNRAS, 332, 945 

Sakamoto T. ct al. , 2005, GCN Circ. 3305 

Sari R., Piran T., 1999, ApJ, 520, 641 

Sari R., Narayan R., Piran T., 1998, ApJ, 497, LI 

Sari R., Piran T., Halpern J., 1999, ApJ, 519, L17 

Stanek K. et al. , 2001, ApJ, 563, 592 

Tagliaferri G. et al. , 2005, Nature, 436, 985 

MacFadyen A., Woosley S., Heger A., 2001, ApJ, 550, 410 

Yost S., Harrison F., Sari R., Frail D., 2003, ApJ, 597, 459 

Zeh A., Klose S., Kann D., 2005, ApJ, accepted (astro- 

ph/0509299) 
Zhang B., Meszaros P., 2001, ApJ, 552, L35 
Zhang B., Kobayashi S., Meszaros P., 2003, ApJ, 595, 950 
Zhang B. et al. , 2005, ApJ, submitted (astro-ph/0508321) 



© 2005 RAS, MNRAS 000, 000-000 



